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Abstract 

The dependence of the Maki-Thompson (MT) and of the density of states de- 
pletion (DOS) contributions from superconducting fluctuations (SF) to NMR- 
NQR relaxation is derived in the framework of the diagrammatic theory, ap- 
plied to layered three dimensional (3D) high T c superconductors. The regular- 
ization procedure devised for the conductivity (Buzdin and Varlamov, Phys. 
Rev. B, 58, 14195 (1998)) is used in order to avoid the divergence of the 
DOS term. The theoretical results are discussed in the light of NMR-NQR 
measurements in YBCO and compared with the recent theory (Eschrig et al., 
Phys. Rev. B 59, 12095 (1999)), based on the assumption of a purely 2D 

spectrum of fluctuations. 
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I. INTRODUCTION 



The normal state of high temperature superconductors (HT C SC) is characterized by un- 
usual properties, most of them still lacking of a comprehensive theoretical description. In 
particular, the transfer of spin excitations from low to high-energy range (spin-gap open- 
ing) well above T c , detected in the generalized susceptibility from Ti and neutron scattering 
measurements (for a recent review see Ref. 1), and the related quasi-particle gap, observed 
in ARPES 2 ' 3 , have been tentatively related to superconducting fluctuations (SF) of various 
nature. Among them, one could mention preformed pairs without long-range phase coher- 
ence, possibly along stripes 4,5 , spin and/or charge density waves 6 ' 7 , coupling of a d-wave 
symmetry order parameter with spin fluctuations 8 ' 9 , order parameter fluctuations well be- 
yond the perturbative approach 10 " 12 and quantum critical point fluctuations 13 . For a review 
on precursor pairing correlations and a survey of the various scenarios see Ref. 14. Further- 
more the magnetic field has been argued 15 to have a role on SF of overdoped HT C SC also, 
by inducing a spin-gap from T c (0) to T C (H). 

The role of the magnetic field is crucial in NMR-NQR attempts to study SF in the 
vicinity of T c + . Indeed, the most direct contribution to SF, namely the Aslamazov-Larkin 
term, responsible of paraconductivity 12 , is not effective in causing an extra-contribution 
to NMR-NQR relaxation. In principle, the SF contribution to the relaxation rates, for 
T — > T c + , are the Maki-Thompson (MT) one (related to the pairing of a carrier with itself 
at a previous stage of motion), and the reduction due to the depletion in the single-particle 
density of states (DOS), when fluctuating pairs are created 12 . These two terms might have a 
different sensitivity to the presence of a magnetic field, which acts as a pair-breaking factor. 

The first NMR experimetal observation 16 of the role played by SF in HT C SC was based 
on the comparison of 63 Cu relaxation rate W in YBCO in the absence (i.e. NQR) and in the 
presence of a magnetic field of about 6 T. Within about 10 degrees above T c (0), W(NQR) 
was found to decrease upon application of the field by a factor about 5-^15%. A qualitative 
interpretation of the experimental observation was given by assuming that the field reduced 
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the MT term to about 25%, while the more robust DOS term was little affected by the field. 
The Equations used in these estimates 16 were the ones 17 pertaining to a three-dimensional 
(3D) layered spectrum of excitations (with anisotropy parameter r = 2£ 2 (0) /d 2 ~ 0.1, where 
£c(0) is the correlation lenght of the Cooper pair at zero temperature, along the c-axis, and 
d the interlayer distance). The occurence of a pure 2D regime of SF could be ruled out, on 
the basis of the absolute value and of the temperature dependence of the SF contribution 
to 63 Cu. 

The first systematic analysis of the field dependence of NMR relaxation rate in YBCO 
was carried out in 1998 by Mitrovic et al. 18 , by varying the field from 2.1 up to 27.3 T. 
The 63 Cu 7\ was probed 18 through the contribution to the 17 echo dephasing. For field in 
the range 6-^8 T, the results derived in this way 18 were found to coincide with the direct 
measurements of 63 Cu T ± . At high field, W DOS was argued to be strongly reduced by the 
field. These data 18 , as well as the field dependence of 17 0(2,3) Knight shift 19 , have been 
interpreted on the basis of a theory for the DOS contribution due to Eschrig et al. 20 , which 
extended analytical approaches 17 to include short wave-lenght and dynamical fluctuations, 
in the assumption of a 2D regime. 

Recently, Gorny et al. 21 reported precise 63 Cu relaxation measurements in YBCO for 
H = 0, 8.8 and 14.8 T, finding no magnetic field dependence in a wide temperature range. 
A possible dependence of the field effect on the amount of doping could be suspected since, 
at the same time, no field effect had been observed 22 in underdoped YBCO. 

In this paper we derive the magnetic field dependence of the MT and DOS contributions 
to NMR Ti, in the framework of a diagrammatic description, for arbitrary values of the re- 
duced field (3 = 2H/H c2 (0), for a 3D layered spectrum of fluctuations, which should pertain 
to the case of YBCO with low anisotropy parameter. In order to remove the logarithmic 
divergence present in the DOS term, here we use the method devised 23 for transverse conduc- 
tivity, in which regularization requirement, analogous to the ones for the nuclear relaxation 
rate, is present. Furthermore, we briefly discuss the role of the long wave-vector fluctuations 
and of the dynamical fluctuations. Our analytical conclusive expressions are compared with 
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the numerical solutions 20 for the 2D regime and with the experimental measurements carried 
out until now. 



II. FIELD DEPENDENCE OF THE SF CONTRIBUTION NUCLEAR 

RELAXATION 

In the following, we extend the diagrammatic theory for the SF contribution to NMR- 
NQR relaxation rate, to include the effects due to the presence of a magnetic field along the 
c-axis, in HT C SC. 

In the presence of the field the MT and DOS contributions to the relaxation rates 12 ' 17 
must be evaluated by starting from the usual inclusion in the momentum q of the term 
(— 2e/c)A, with A = §(— yi + xj). The integration over q, in the ab plane, is substituted 
by a sum over the Landau levels. Thus the MT and DOS contributions 17 to W assume the 
forms 

W MT lo . it (3 



-(P,e) — - 7 • 

8E f t (e - lv ) V L ^ + (3( n + 1/2) ySv + P(n + 1/2) + r 

1 



+ 



y/e + J{n + l/2)yje + 0(n + 1/2) + r 



(1) 



and 



w dos n i//3 1 

v -l— — (/3, 6 ) = -JL- . k (t t ) ■ V . . (2) 

where W° is the ordinary relaxation rate in the absence of SF (the Korringa one in a Fermi 
gas-like model). In the above Equations, E F is the Fermi energy, r the single particle 
collision time, 7^ = ^q/Dt^, (with D = Epr/m, 2D carrier diffusion constant and Tit~ 1 
depairing energy) is a dimensionless factor which, in the limit B —> 0, takes into account 
the pair-breaking effect and e = (T — T c )/T c is the reduced temperature. k(Tt) in Eq. (2) 
is the function 

K(TrH rc(3, 1 _/*P>, *v<, 



7T 4ttTt [^(1/2) - ip(l/2 + 1/4ttTt)] + ^'(1/2) 
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4Tr, Tr > 1 



where, K B — h — 1 and ip is the Euler digamma function. The expressions (1) and (2) are 
general and hold for any magnetic field, provided that (3 = 2B/B c2 (0) <C 1 (see Ref. 24). 
The DOS term shows a logarithmic divergence in the sum. 

The behaviours of W MT and W DOS , for small fields, can easily be derived from Eqs. (1) 
and (2), by a straightforward expansion. The intermediate and strong field regimes, namely 
e <C (3 <C 1, require some care. For the MT contribution (Eq. (1)) the sum over n converges 
rapidly and thus, for strong fields, only the term n = can be taken into account. For the 
DOS contribution (Eq. (2)), however, the logarithmic divergence must be removed. This 
can be achieved by considering the cut-off independent difference 



AW UUb ((3, e) = W UUb ((3, e) - W UUb (0, e) (3) 
The zero field expression is rewritten in the form 

WD °y>* = ■ 2 K (Tr) lim £ In CT+g+y/FTrWg (4) 

W° E F r V ; M n tJ y/F+]3n + y/e + r + Pn V ; 

and then, Eq. (3) becomes 

AW DOS ((3,s) h ™( + (3n + [3 + y/e + r + fin + [3 



n 2n(Tr) f ln Ve + [Jn + [3 + ye + r + (3n + (3 

g } (5) 

2 yje + (3n + (3/2^e + r + f3n + (3/2> 



where the summation has been extended up to n — > oo, since AW DOS (j3, e) now involves a 
summation which is convergent (the n-th term being proportional to n _3//2 for large n). 

Now we are going to discuss the behaviours of W DOS and W MT in the asymptotic limits, 
corresponding to different temperature and field regimes. 



A. Weak magnetic fields 

In the weak field regime, namely (3 <ti e, Eqs. (1), (2) can be expanded in powers of (3. 
For the DOS term, in particular, the Euler-MacLaurin formula 

E /(") = [ N f(x)dx + kf(N) + /(0)] + hf'(N) - /'(0)] + . . . 



yields 
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B. Intermediate and strong field regimes 

For the MT term in the strong field regime (f3 3> max{e, r, 7^}), by expanding Eq. (1) 
in powers of one has 



TMT 



■(e,7 v ,r < /3) 



3vr 3 1 



(8) 



W° v ~' " r ' IQE f t{3 

In the intermediate case (namely e, 7^ <C /3 <C r and e <C /5 <C 7^, r), which can become 
relevant for a 3D layered compound, the series expansions in terms of the smallest parameters 
yield 



W MT IT 1 

-(e,7 v ,«/?«r) = 4.57- 



16E F r 



(9) 



^yMT 7T 1 , v /max{ 7 ^,r} 
(£</?« 7^,r) = — In v 



(10) 



For the field dependence of the DOS correction to the zero field contribution (Eq. (5)), 
in strong field regime, one can take into account the n = term only: 



A^° 5 (/i»max{e,r}) _ h 2/<Tr) j ln . ^ 



(11) 



W° EpT ' / V"e(^E + v^+7)J 

The n > 1 terms, neglected in this evaluation, yield a correction of 0.02 in the bracket term 
in Eq (11). 



6 



In the intermediate fields, namely for e <C (3 <C r (corresponding to a 3D layered regime 
of fluctuations), by means of an expansion of Eq. (5) in terms of /3/2r one has 

W° E f t v ' V 2r v ; 



III. DISCUSSION AND COMPARISON WITH EXPERIMENTAL RESULTS 

As it appears from the theoretical treatment given in the previous Section, the depen- 
dence on the magnetic field of the SF contributions to the NMR-NQR relaxation rate is 
a rather delicate issue, because of the non-trivial interplay of several parameters, such as 
(T — T c ), the reduced field f3 = 2B/B c2 (0), the anisotropy parameter r, the elastic collision 
time r and the anelastic phase-breaking time t v . 

Eschrig et a/. 20 , in their numerical extension of the previous analytical approaches 17 , have 
considered arbitrary values of k(Tt) and taken into account short wave-lenght and dynamical 
fluctuations. However, this generalization has required the restriction to a purely 2D regime 
of SF, which seems questionable in view of experimental findings in YBCO, pointing out a 
crossover to 3D fluctuations well above T c , at least for relatively small fields 25-27 . On the 
other hand, it could be remarked that dynamical fluctuations are relevant only when the 
field is comparable to B c2 (0) ~ 100 4- 120 T. Therefore the reduction to static and small 
wave-vector fluctuations should not invalidate our conclusive expressions given above. 

To illustrate the theoretical expressions, derived in Section II, we plot in Figs. 1 and 
2 the temperature and field behaviours of W DOS and W MT , with a choice of parameters 
appropriate to YBCO optimally doped. 

In Fig. 3 the experimental results from various authors (and from different YBCO 
samples, about optimally doped) are compared with the behaviour expected for DOS con- 
tribution, according to Eqs. (3) and (5). A relatively small field dependence, of both DOS 
and MT contributions, up to a reduced field of about 0.2, is noticed. For strong fields a 
reduction of the MT term can be expected, while the DOS term seems only slightly affected 
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(see Figs, lb and 2b). As a consequence, either for s-symmetry of the fluctuating Cooper 
pair (namely, presence of both the MT and DOS effects of SF on Ti) or for d-symmetry 
(and therefore, no MT contribution) only a slight dependence of Ti, on the magnetic field, 
should be detected, as it is illustrated in Fig. 3, for the DOS term. In particular the exper- 
imental data 18 at strong fields ((3 > 0.2) can hardly be justified. The discrepancy between 
the theoretical description and the experimental findings, at strong fields, might be related 
to the general framework involved in our treatment. Howerver, it should be stressed that 
a breakdown of the Fermi liquid picture should not invalidate the expression for the DOS 
contribution, which is independent from the pairing mechanism and the normal state prop- 
erties. The dimensional crossover (3D— >2D) for strong fields is taken into account in our 
equations, but it could be argued that dynamical and short wave-lenght fluctuations are no 
longer negligible in a 2D strong field regime. Finally, it should be taken into consideration 
the possibility that the method of an indirect estimate of 63 Cu Ti from 17 echo dephasing 18 
could be invalidated in strong fields. 

Summarizing, from the behaviour of 63 Cu 

T 16,18,19,21 

one can infer that a DOS SF con- 
tribution to the nuclear relaxation is present in the vicinity of T c . The MT term is more 
elusive, requiring a s-wave component in the spectrum of the SF and being possibly strongly 
sample-dependent, through the pair-breaking and impurities effects. The field dependence 
of DOS term is a rather delicate issue, complicated by various crossovers and parameters. 
Further theoretical and experimental work is required for a firm conclusion. 
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CAPTIONS FOR FIGURES 



Fig. 1: Theoretical behaviours for the DOS contribution to the nuclear spin-lattice re- 
laxation rate (normalized to the value in the absence of SF, W°), according to Eqs. (3) and 
(5), in the text, as a function of temperature for different values of the field (a), and as a 
function of the reduced magnetic field (3 = 2B/B c2 (0) for fixed values of the temperature 
(b). The curves have been derived in correspondence to choices of the upper critical field 
B C 2(0) ~ 120 T, the Fermi energy E F = 3500 K B and the single-particle collision time 
t = 2 • 10- 14 s. 

Fig. 2: Theoretical behaviours for the MT contribution to the nuclear spin-lattice relax- 
ation rate, according to Eq. (1), in the case of strong pair-breaking (7^ ~ 0.3). Part (a) of 
Fig. 2 shows W MT ((3,e)/W° as a function of temperature (for different values of the field) 
while part (b) reports it as a function of the reduced field, (5 = 2B / B c2 (0), in correspondence 
to fixed values of the temperature. The choice of parameters is the same as in Fig. 1. 

Fig. 3: Comparison between the experimental results, obtained at T ~ 95 K (□ Mitrovic 
et al. 18 ; O Gorny et al. 21 ; □ da Carretta et al. 16 ; x da Carretta et al. (data unpublished)) 
and the theoretical predictions for the DOS term as a function of the reduced field in 
correspondence to the usual choice of parameters (Figs. 1 and 2). The experimental data 
by Gorny et al. 22 (O) and by Carretta et al. (x) (unpublished) have been reported in 
correspondence to the value W DOS (0,e) ~ — 0.16W^°, in order to analyze their possible field 
dependence. 
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